Measurement of the tt Production Cross Section in pp Collisions at a/s = 1-96 TeV 
using Missing Ex+jets Events with Secondary Vertex 6- Tagging 
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We present a measurement of the tt production cross section in pp collisions at ^/s — 1.96 TeV which 
uses events with an inclusive signature of significant missing transverse energy and jets. This is the 
first measurement which makes no explicit lepton identification requirements, so that sensitivity 
to W ^ Tiy decays is maintained. Heavy fiavor jets from top quark decay are identified with a 
secondary vertex tagging algorithm. From 311 pb~^ of data collected by the Collider Detector at 
Fermilab we measure a production cross section of 5.8 db 1.2(stat.)lQ"7(syst.) pb for a top quark 
mass of 178 GeV/c^, in agreement with previous determinations and standard model predictions. 

PACS numbers: 13.85.Ni, 14.65. Ha 



At the Tevatron pp collider top quarks are produced 
mainly in pairs through quark- ant iquark annihilation and 
gluon-gluon fusion processes. In the standard model 
(SM) the calculated cross section for pair production is 
Q.ltoi pb JL] for a top mass of 178 GeV/c^ 2], and 
varies linearly with a slope of —0.2 pb/GeV/c^ with the 



top quark mass in the range 170 GeV/c^ < < 190 
GeV/c^. Because the Cabibbo-Kobayashi-Maskawa ma- 
trix element Vtb is close to unity and rrit is large, the SM 
top quark decays to a W boson and a b quark almost 
100% of the time. The final state of top quark pair pro- 
duction thus includes two W bosons and two 6-quark jets. 
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When only one W decays leptonically, the tt event typi- 
cahy contains a charged lepton, missing transverse energy 
from the undetected neutrino, and four high transverse 
energy jets, two of which originate from b quarks Pre- 
vious cross section analyses 4, 5, 6] select this distinct tt 
signature by requiring well-identified leptons (e, ja) with 
high transverse momentum. In this Letter we describe a 
tt production cross section measurement which is sensi- 
tive to leptonic W decays regardless of the lepton type, 
and has a sizable acceptance to r decays of the W bo- 
son. The direct identification of r from W ^ tv decays 
suffers from a very low efficiency, thus our measurement, 
using data collected by a multijet trigger, selects top de- 
cays by inclusively requiring a high-p^ neutrino signature 
rather than charged lepton identification. Events with 
well-identified high-p^ electrons or muons are explicitly 
vetoed in order to avoid statistical overlap and provide 
complementary results with respect to lepton-based mea- 
surements. 

Results reported in this Letter are obtained from 311 
pb~^ of data from pp collisions at ^/s = 1.96 TeV 
recorded by the Collider Detector at Fermilab (CDF II). 
The CDF II detector has been described in detail else- 
where 0. It consists of a magnetic spectrometer sur- 
rounded by calorimeter systems and muon chambers. 
The momenta of charged particles are measured up to a 
pseudorapidity of |7^| = 1.0 in the central tracking cham- 
ber, which is inside a 1.4 T superconducting solenoidal 
magnet. Microstrip silicon vertex detectors, located im- 
mediately outside the beampipe, provide precise track 
reconstruction useful for vertexing and extend the |7^| 
coverage of the tracking system up to |7^| = 2.0. Electro- 
magnetic and hadronic sampling calorimeters, arranged 
in a projective- tower geometry, surround the tracking 
systems and measure the energy and direction of elec- 
trons, photons, and jets, providing good hermeticity in 
the pseudorapidity range \r]\ < 3.6. In addition, the 
calorimeter system allows the detection of high-p^ neutri- 
nos by the measurement of the missing transverse energy, 
3] . Muon systems reside outside the calorimeters and 
consist of layers of drift chambers which allow the re- 
construction of track segments for penetrating particles. 
The beam luminosity is determined using gas Cherenkov 
counters (CLC) located in the region 3.7 < \r]\ < 4.7 
which measure the average number of inelastic pp colli- 
sions per bunch crossing. The uncertainty on the lumi- 
nosity measurement, due to CLC acceptance and theoret- 
ical uncertainties on the total inelastic pp cross section, 
is 6% 0. 

The data sample used in this analysis is collected by a 
multijet trigger which requires four or more Et > 15 
GeV clusters of contiguous calorimeter towers, and a 
total transverse energy clustered in the calorimeter of 
^ Et > 125 GeV. The initial data sample consists of 4.2 
million events. The understanding of signal acceptances 
and efficiencies relies on a detailed simulation of the pro- 



duction processes and the detector response. Inclusively 
decaying tt events, assuming a top quark mass of 178 
GeV/c^, are simulated using pythia v6.2 and her- 
WIG v6.5 10] generators in conjunction with the cteq5l 
[Til parton distribution functions, QQvQ.l 
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for the 

modeling of b and c hadron decays, and a full simula- 
tion of the CDF II detector. Jets are identified as groups 
of calorimeter tower energy deposits which fall within a 
cone of radius AR = y^^^^^^ ^ 0.4. Jet energies, 
after the absolute energy scale setting, are corrected for 
calorimeter non-linearity, losses in the gap between tow- 
ers, and multiple interactions |l^ . 

The tt signature used in the present study (^^+jets) 
consists of large missing transverse energy, associ- 
ated with the neutrino from the leptonic decay of a W bo- 
son, and jets. Since the resolution, cr(^^), is observed 
to degrade as a function of the total transverse energy of 
the event, in the form cr(^^) ex Et ^14], the miss- 
ing Et significance, defined by = \^t\I ^t-, is 
used for event selection. Besides the requirement of large 
further background rejection can be achieved by ex- 
ploiting ^^-related geometrical properties. The neutrino 
direction for tt ^t + J^^s events is in general uncor- 
related with the jet directions in the transverse plane. 
On the contrary, background events for which the $t is 
mainly due to jet energy mis- measurement or to 6-quark 
semileptonic decays are more likely to have I$t aligned 
with a jet direction. The minimum Acj) difference be- 
tween the ^T j^^ i^ event, min A(j){^rj.^ jets)^ 
is used as an analysis cut. Events containing identi- 
fied high-pT electrons or muons, as defined in 5], are 
removed, in order to increase the relative contribution 
oi W ^ TV decays and to provide a statistically inde- 
pendent sample with respect to other lepton+jets cross 
section analyses. 

To reject events with only light quark or gluon jets, 
we require at least one jet to be identified as originat- 
ing from a 6-quark. The presence of 6-jets ("tags") 
is established by the identification of secondary decay 
vertices using the secvtx algorithm 5]. An optimiza- 
tion procedure, aimed at minimizing the relative statis- 
tical uncertainty on the cross section measurement, sets 
Njet{ET > 15 GeV; \r]\ < 2.0) > 4, ^^'^ > 4.0 GeV 
and min A(j){^ J. ^ jets) > 0.4 rad as the best kinematical 
selection cuts. After these requirements the data sample 
is reduced to 597 events with an expected pre-tag signal 
to background ratio S/B ~ 1/5. 

Background events with 6-tags arise from QCD heavy 
flavor production, electroweak production of W bosons 
associated with heavy flavor jets, and from false identi- 
flcation by the secvtx algorithm. The overall number 
of background 6-tags is estimated from the multijet sam- 
ple by applying a parametrization of the per-jet tagging 
probability. The tagging probability is calculated using 
~ 879, 000 multijet data events with exactly three jets 
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having Et > 15 GeV and |?7| < 2.0. It is parametrized 
on a per-jet basis as a function of the jet Et^ track mul- 
tiphcity, and the projection of the along the jet di- 
rection, defined by = • cos A0(^^, jet). The tt 
contamination in this control sample is negligible. The 
jet 6-tagging rate, calculated as the ratio of the number of 
6-tagged jets to the number of jets with at least two good 
tracks with hits in the silicon detector, is shown in Fig.^ 
The parametrization accounts for background 6-tags 
from bb production processes, whose 6-tag rate is en- 
hanced at high positive values of the due to the cor- 
relation of and 6-jet direction in the case of a semilep- 
tonic 6-quark decay. The extrapolation of the 3-jet 6-tag 
rate to higher jet multiplicity events is checked by com- 
paring the predicted and observed 6-tag rates as a func- 
tion of the number of jets in the multijet sample without 
the kinematical selection on J^^^ and min A^(^^, jets), 
as shown in Fig.I2fa). The capability of the parametriza- 
tion to track possible sample composition changes in- 
troduced by the kinematical selection is established by 
constructing the ratio of observed over expected 6-tags 
in two separate control samples of multijet data, de- 
pleted of signal contaminations, as shown in Fig.EJb) and 
Fig.[2fc), respectively. The first control sample is defined 
by ^^'^ < 3.0 GeV^/^ min A0(^^, jets) > 0.3 rad; the 
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FIG. 1: Tagging rate probability parametrization in the mul- 
tijet sample as a function of jet Et (a), track multiplicity (b), 
and missing Et projection (c). 
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FIG. 2: Ratio of observed to expected 6-tags as a function of 
the number of jets in the data before the optimized kinemati- 
cal selection on and min A(l){^rp, jets) (a), and in control 
samples: < 3.0 GeV^/^ min A0(^^, jets) > 0.3 rad (b); 
and > 3.0 GeV^/^ min A0(^^, jets) < 0.3 rad (c). 



second as ^^'^ > 3.0 GeV^/^ min A^(^^, jets) < 0.3 
rad. The 6-tag rate parametrization is found to predict 
the non-tt background to within 10% in the 4 < Njet < 6 
region where 96.4% of the tt signal is expected. This 
value is thus adopted as the systematic uncertainty on 
the background prediction. Additional checks, using low 
jet multiplicity events from leptonic W data samples, 
show agreement between observed and predicted 6-tags 
to within the quoted uncertainty. 

The sample selected with the optimized kinematical 
requirements described above and the requirement of at 
least one 6-tagged jet consists of 106 events, containing 
a total of 127 6-tagged jets. The number of 6-tagged 
jets yielded by background processes in that sample is 
expected to be Nexp = 67.4 ± 2.7 ± 6.7. The first uncer- 
tainty contribution is due to the limited number of events 
in the 3-jet sample used for the tagging rate parametriza- 
tion, while the second contribution is the 10% system- 
atic uncertainty on the 6-tag rate parameterization dis- 
cussed above (Fig. [2j. The 597 pre-tag data events are 
expected to contain a non- negligible tt component, due to 
the S/N enhancement provided by the kinematical selec- 
tion, which yields a background overestimate. The N^xp 
value is corrected, to N'^^^ = 57 A ±8.1, to account for 
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TABLE I: Number of 6-tagged jets expected from Monte Carlo ti production using a{tt) = 6.1 pb {rrit = 178 GeV/c^), 
predicted by the tagging rate parametrization, and observed in the selected sample, as a function of the jet multiplicity. The 
total uncertainty on the number of predicted background 6-tags is the sum in quadrature of the statistical uncertainty and of 
a 10% systematic uncertainty. The number of background 6-tags corrected for the tt contamination in the pre-tagging data 



sample is also provided for the signal re^ 


^ion. Uncertainties on si^ 


^nal contributions 


are statistical only. 




Number of jets 


3 


4 


5 


> 6 


tt ee 


0.08 ±0.01 


0.41 ±0.03 


0.18 ±0.02 


0.04 ±0.01 


tt e/i 


0.06 ±0.01 


0.29 ±0.02 


0.11 ±0.01 


0.05 ±0.01 


tt iifi 


0.01 ±0.01 


0.05 ±0.01 


0.01 ±0.01 


0.01 ±0.01 


tt er 


0.11 ±0.01 


0.93 ± 0.04 


0.38 ± 0.03 


0.15 ±0.02 


tt IIT 


0.05 ±0.01 


0.29 ±0.02 


0.15 ±0.02 


0.06 ±0.01 


tt TT 


0.06 ±0.01 


0.58 ± 0.03 


0.26 ±0.02 


0.05 ±0.01 


ti ^ e-\- jets 


0.68 ± 0.04 


6.61 ±0.11 


8.70 ±0.13 


4.25 ±0.09 


tt ^ /J -\- jets 


1.07 ±0.04 


11.92 ±0.15 


6.56 ±0.11 


2.47 ±0.07 


ti ^ T -\- jets 


1.00 ±0.04 


10.98 ±0.14 


11.71 ±0.15 


5.53 ±0.10 


ti jets 


0.01 ±0.01 


0.09 ±0.01 


0.14 ±0.02 


0.22 ±0.02 


ti^ X 


3.13 ±0.08 


32.15 ±0.24 


28.14 ±0.23 


12.83 ±0.15 


Background 6-tagged jets 


32.68 ± 3.46 


37.53 ±4.14 


21.44 ±2.76 


8.47 ± 1.40 


Corrected background 6- tagged jets 




33.14 ±4.01 


17.58 ±2.85 


6.71 ±2.78 


Observed 6- tagged jets 


31 


53 


55 
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In order to further establish the ti signal in the se- 
lected data we perform binned likelihood fits to kinemat- 
ical distributions. The stability of the fitting technique 
is checked using simulations with known signal fractions. 
We fit and A(j){l^j.^ tagged jet) data distributions to 
the sum of a ti and a background template. The former 
is obtained from Monte Carlo ti inclusive events; the lat- 
ter is derived from the tagging rate parametrization ap- 
plied to the data. Results for the A(j){l^rj.^ tagged jet) 
fit are displayed in Fig. [31 The fitted ti components in 
the selected data (68 ± 12% and 44 ± 12% for the 
and A(j){^j.^ tagged jet) fits respectively) are in agree- 
ment with the overall prediction calculated from 6-tag 
counting, before any correction to account for the ti con- 
tamination in the pre-tagging data sample (47±5% deter- 
mined comparing the number of expected and observed 
6-tags, for Njet > 4, in Table HI). 

The efficiency of the trigger, the kinematical selec- 
tion, and the 6-tagging algorithm, are evaluated using 
inclusive ti Monte Carlo events. The combined effi- 
ciency of the trigger and kinematical selection amounts 
to Ckin = 4.88% ± 0.43% for a top mass of rrit = 178 
GeV/c^, where the dominant uncertainty is determined 
by the comparison of the results from the PYTHIA and 
HERWIG generators. Other sources of systematic uncer- 
tainty are evaluated by varying Monte Carlo generation 
settings with respect to the default values, and by ap- 
plying systematic variations of the jet correction factors; 
the trigger acceptance uncertainty is determined by com- 
paring trigger turn-on curves between Monte Carlo and 
data events (Tabled). The average number of 6-tags per 
Monte Carlo ti event is found to be e^^^ = 0.789 ± 0.046, 
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FIG. 3: A(l){^r^, tagged jet) distribution for data after kine- 
matical selection and with at least one ^-tagged jet. The data 
are fit to the sum of ti signal and background templates as 
described in the text. 



this effect. Table U summarizes the number of 6-tagged 
jets expected from Monte Carlo simulation for each ti 
decay mode satisfying the kinematical and > 1 6-tag re- 
quirements, as well as the predicted and observed 6-tags 
in the selected data sample as a function of the jet mul- 
tiplicity of the events. In Tabled the numbers of 6- tags 
from events with only three jets are provided as a cross- 
check of the background. The excess in the number of 
6-tags, for N{jet) > 4, is ascribed to top pair production. 
Final states with W ^ rv decays account for ~ 44% of 
the signal acceptance. 
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TABLE II: Relevant sources of systematic uncertainty. 



Source relative error 

Cfcin systematics 

Generator dependence 8.2 % 

Trigger acceptance 2.0 % 

ISR/FSR 2.0 % 

PDFs 1.6 % 

Jet Energy Scale 1.5 % 
Others 

Background prediction method 10.0 % 

Luminosity measurement 6.0 % 

Ctag (SECVTX scale factor) 5.8 % 
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and it is corrected according to the Monte Carlo secvtx 
efficiency scale factor of 0.909 ± 0.060 to reproduce the 
data 6-tagging efficiency 0. 

The cross section is calculated with a Poisson likeli- 
hood function in which the maximum likelihood solu- 
tion for (7(tt) is given by: (j(tt) = — "^""a^e 7^ ? where 

Nobs and A^exp number of 6-tagged jets observed 

and expected from background events by the tagging 
rate parametrization in the selected data, and Cint is 
the integrated luminosity of the multijet data sample. 
The input parameters: Cint, ^kim ^tag^ ^exp 
subject to Gaussian constraints. With these input val- 
ues we measure a top pair production cross section of 
5.8 ± 1.2(stat.)lo.7(sys^-) Pt>- Additional samples of in- 
clusive tt Monte Carlo events generated with different mt 
values in the range [130,230] GeV/c^ are used to com- 
pute the cross section measurement dependence on m^. 
The cross section measurement changes by ±0.05 pb for 
each ^1 GeV/c^ change in the top quark mass from the 
initial value of 178 GeV/c^. For instance, we measure 
a{tt) = 6.0 ± 1.2(stat.)lQ7(syst.) pb assuming a top 
quark mass of 175 GeV/c^. The change is due to the 
varying signal selection efficiency with top quark mass. 

In conclusion we report the first measurement of the 
top pair production cross section of a{tt) = 5.8lJ;4 pb 
using an inclusive selection of ^2^+ jets tt decays. The 
result is complementary to other cross section analyses 13, 
1^ , maintains high sensitivity with respect to W ^ ru 
tt decays, and is in good agreement with SM calculations 
and previous measurements. 
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